- We load the CITE-seq human PBMC data. You can follow the Seurat tutorial or download the analyzed R object.

- We run WLR on the RNA data.
ptm <- proc.time()
cbmc_rna <- RunWLR(object = cbmc_rna, n_iters = 1000, wlr_pca_dims = 20, wlr_k = 30, ratio = 0.1, norm_type = "l1")
Run Bootstrapping
Finalize
ptm1 <- proc.time()
print(x = ptm1 - ptm)
user system elapsed
42.979 13.618 70.487
- We run WLR on the ADT data.
print(x = ptm1 - ptm)
user system elapsed
34.798 11.520 58.947
- We fuse the RNA and ADT data. We see that CD4+ and CD8+ T cells are more defined by the ADT data in the surface protein space, while Dendritic cells (DCs) are more defined in the RNA-seq data due to the absence of its surface marker in the ADT data.

- We perform UMAP dimensional reduction on the fused data. Additionally, clustering can be performed as well.

LS0tCnRpdGxlOiAiV0xSIFR1dG9yaWFsIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCjEuIFdlIGxvYWQgdGhlIFtDSVRFLXNlcV0oaHR0cHM6Ly93d3cubmF0dXJlLmNvbS9hcnRpY2xlcy9ubWV0aC40MzgwKSBodW1hbiBQQk1DIGRhdGEuIFlvdSBjYW4gZm9sbG93IHRoZSBbU2V1cmF0IHR1dG9yaWFsXShodHRwczovL3NhdGlqYWxhYi5vcmcvc2V1cmF0L3YzLjEvbXVsdGltb2RhbF92aWduZXR0ZS5odG1sKSBvciBkb3dubG9hZCB0aGUgYW5hbHl6ZWQgW1Igb2JqZWN0XShodHRwczovL2RyaXZlLmdvb2dsZS5jb20vZHJpdmUvZm9sZGVycy8xaG12ZjF2SmRHMkRONjBmbDFqekxXbUNNLVNOd1lqalg/dXNwPXNoYXJpbmcpLgoKYGBge3Igc2V0dXAsIGZpZy5oZWlnaHQ9NSwgZmlnLndpZHRoPTEwfQpsaWJyYXJ5KFdMUikKbGlicmFyeShTZXVyYXQpCmNibWMgPC0gcmVhZFJEUyhmaWxlID0gImNibWMuZGF0YS5yZHMiKQpDb21iaW5lUGxvdHMocGxvdHMgPSBsaXN0KERpbVBsb3QoY2JtYywgcmVkdWN0aW9uID0gInVtYXAiLCBncm91cC5ieSA9ICJybmFDbHVzdGVySUQiLCBsYWJlbCA9IFRSVUUpICsgZ2dwbG90Mjo6Z2d0aXRsZSgiUk5BLXNlcSBEYXRhIiksIERpbVBsb3QoY2JtYzEsIHJlZHVjdGlvbiA9ICJ0c25lX2FkdCIsIGdyb3VwLmJ5ID0gInJuYUNsdXN0ZXJJRCIsIGxhYmVsID0gVFJVRSkgKyBnZ3Bsb3QyOjpnZ3RpdGxlKCJBRFQgU3VyZmFjZSBQcm90ZWluIERhdGEiKSkpCmBgYAoKMi4gV2UgcnVuIFdMUiBvbiB0aGUgUk5BIGRhdGEuCmBgYHtyIFJOQX0KZGF0YV9ybmEgPC0gY2JtY0Bhc3NheXMkUk5BQGNvdW50c1tjYm1jQGFzc2F5cyRSTkFAdmFyLmZlYXR1cmVzLF0KY2JtY19ybmEgPC0gSW5pdGlhdGVXTFIocmF3X2RhdGEgPSBkYXRhX3JuYSwgcHJvamVjdF9uYW1lID0gIlJOQSIsIG5vcm1hbGl6YXRpb24gPSAiY29zaW5lIiwgcGNhX2RpbXMgPSAyMCwga2VybmVsID0gImdhdXNzaWFuIiwgayA9IDMwKQpjYm1jX3JuYSA8LSBHZW9tU2tldGNoKG9iamVjdCA9IGNibWNfcm5hLCBnZW9tX3BjYV9kaW1zID0gMjAsIGdlb21fc2l6ZSA9IDEwMDAsIHNrZXRjaF9uX3BjYSA9IDIwLCBza2V0Y2hfayA9IDMwKQpwdG0gPC0gcHJvYy50aW1lKCkKY2JtY19ybmEgPC0gUnVuV0xSKG9iamVjdCA9IGNibWNfcm5hLCBuX2l0ZXJzID0gMTAwMCwgd2xyX3BjYV9kaW1zID0gMjAsIHdscl9rID0gMzAsIHJhdGlvID0gMC4xLCBub3JtX3R5cGUgPSAibDEiKQpwdG0xIDwtIHByb2MudGltZSgpIApwcmludCh4ID0gcHRtMSAtIHB0bSkKYGBgCgozLiBXZSBydW4gV0xSIG9uIHRoZSBBRFQgZGF0YS4KYGBge3IgQURUfQpkYXRhX3BydCA8LSBjYm1jQGFzc2F5cyRBRFRAY291bnRzCmNibWNfcHJ0IDwtIEluaXRpYXRlV0xSKHJhd19kYXRhID0gZGF0YV9wcnQsIHByb2plY3RfbmFtZSA9ICJQcm90ZWluIiwgbm9ybWFsaXphdGlvbiA9ICJjIiwga2VybmVsID0gImdhdXNzaWFuIiwgayA9IDMwLCB2ZXJib3NlID0gRkFMU0UpCmNibWNfcHJ0IDwtIEdlb21Ta2V0Y2gob2JqZWN0ID0gY2JtY19wcnQsIGdlb21fc2l6ZSA9IDEwMDAsIGdlb21fcGNhX2RpbXMgPSAzLCBza2V0Y2hfbl9wY2EgPSAwLCBza2V0Y2hfayA9IDMwKQpwdG0gPC0gcHJvYy50aW1lKCkKY2JtY19wcnQgPC0gUnVuV0xSKG9iamVjdCA9IGNibWNfcHJ0LCBuX2l0ZXJzID0gMTAwMCwgcmF0aW8gPSAwLjEsIHdscl9rID0gMzAsIG5vcm1fdHlwZSA9ICJsMSIpCnB0bTEgPC0gcHJvYy50aW1lKCkgCnByaW50KHggPSBwdG0xIC0gcHRtKQpgYGAKCjQuIFdlIGZ1c2UgdGhlIFJOQSBhbmQgQURUIGRhdGEuIFdlIHNlZSB0aGF0IENENCsgYW5kIENEOCsgVCBjZWxscyBhcmUgbW9yZSBkZWZpbmVkIGJ5IHRoZSBBRFQgZGF0YSBpbiB0aGUgc3VyZmFjZSBwcm90ZWluIHNwYWNlLCB3aGlsZSBEZW5kcml0aWMgY2VsbHMgKERDcykgYXJlIG1vcmUgZGVmaW5lZCBpbiB0aGUgUk5BLXNlcSBkYXRhIGR1ZSB0byB0aGUgYWJzZW5jZSBvZiBpdHMgc3VyZmFjZSBtYXJrZXIgaW4gdGhlIEFEVCBkYXRhLgpgYGB7ciBmdXNlLCBmaWcuaGVpZ2h0PTUsIGZpZy53aWR0aD0xMH0KZnVzZWQgPC0gRnVzZVdMUihjYm1jX3JuYSwgY2JtY19wcnQsIHByb2plY3RfayA9IDEwKQpjYm1jQG1ldGEuZGF0YVsiUk5BX1Njb3JlIl0gPC0gZnVzZWQkZnVzZWRfd2VpZ2h0WyJSTkEiLF0KY2JtY0BtZXRhLmRhdGFbIlByb3RlaW5fU2NvcmUiXSA8LSBmdXNlZCRmdXNlZF93ZWlnaHRbIlByb3RlaW4iLF0KQ29tYmluZVBsb3RzKHBsb3RzID0gbGlzdChWbG5QbG90KG9iamVjdCA9IGNibWMxLCBmZWF0dXJlcyA9ICJSTkFfU2NvcmUiLCBzb3J0ID0gVFJVRSksIFZsblBsb3Qob2JqZWN0ID0gY2JtYzEsIGZlYXR1cmVzID0gIlByb3RlaW5fU2NvcmUiLCBzb3J0ID0gVFJVRSkpKQpgYGAKCjUuIFdlIHBlcmZvcm0gVU1BUCBkaW1lbnNpb25hbCByZWR1Y3Rpb24gb24gdGhlIGZ1c2VkIGRhdGEuIEFkZGl0aW9uYWxseSwgY2x1c3RlcmluZyBjYW4gYmUgcGVyZm9ybWVkIGFzIHdlbGwuCmBgYHtyIHVtYXAsIGZpZy5oZWlnaHQ9NCwgZmlnLndpZHRoPTV9CnVtYXAuZGltcyA8LSB1d290Ojp1bWFwKFggPSBhcy5tYXRyaXgoZnVzZWQkZnVzZWRfZGlzdCksIG5fY29tcG9uZW50cyA9IDIpCmRpbW5hbWVzKHggPSB1bWFwLmRpbXMpIDwtIGRpbW5hbWVzKHggPSBjYm1jQHJlZHVjdGlvbnMkdW1hcEBjZWxsLmVtYmVkZGluZ3MpCmNibWNAcmVkdWN0aW9ucyR1bWFwQGNlbGwuZW1iZWRkaW5ncyA9IHVtYXAuZGltcwpEaW1QbG90KG9iamVjdCA9IGNibWMsIGdyb3VwLmJ5ID0gInJuYUNsdXN0ZXJJRCIsIHJlZHVjdGlvbiA9ICJ1bWFwIiwgbGFiZWwgPSBUUlVFKQpgYGAK